Low Temperature Properties of the Random Field Potts Chain 
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The random field (/-States Potts model is investigated using exact groundstates and finite- 
temperature transfer matrix calculations. It is found that the domain structure and the Zeeman 
energy of the domains resembles for general q the random field Ising case (q = 2). This is also the 
expected outcome based on a random- walk picture of the groundstate. The domain size distribution 
is exponential, and the scaling of the average domain size with the disorder strength is similar for 

q arbitrary. The zero-temperature properties are compared to the equilibrium spin states at small 
temperatures, to investigate the effect of local random field fluctuations that imply locally degenerate 
regions. The response to field perturbations ('chaos') and the susceptibility are investigated. In 
particular for the chaos exponent it is found to be I for q — 2, . . . ,5. Finally for q — 2 (Ising case) 
the domain length distribution is studied for correlated random fields. 

PACS numbers: 05.40.-a, 05.50.+q, 75.50.Lk 



I. INTRODUCTION 

Adding disorder to random magnets creates new ef- 
fects one the prime examples being the random field Ising 
model (RFIM) §, |. Since the randomness couples di- 
rectly to the order parameter one has interesting physics 
in the groundstate (GS) properties || |], [| |), already. 
With a finite temperature the disordering effects of en- 
tropy and local fields compete. 

The purpose of this article is to explore the gener- 
alization or extension of the one-dimensional RFIM to 
the one-dimensional random-field Potts model. The Id 
RFPM is defined on a chain of L spins (with periodic 
boundary conditions being used here) where each spin 
site can be in one of q states. The Hamiltonian is given 
by 



H 



(i) 



where are nearest-neighbour pairs and <5<j !Cr < the 

Kronecker-delta, i.e., 5 CT>cr ' = 1 for a = a' and b G ^i = 
otherwise . The vector on is a unit vector pointing ran- 
domly in one of the q spin directions and hi is the local 
field strength which is chosen either constant or randomly 
distributed. In the numerical computations presented be- 
low we have used a Gaussian distribution for the random 
fields with (hi) = and (hi 2 ) — h r 2 , such that h r /J 
is a measure for the strength of the disorder, except at 
finite T where fields are uniformly distributed between 
[—<5, +5}. In the following we denote the one-dimensional 
versions of the RFIM and the RFPM with RFIC (random 
field Ising chain) and RFPC(random field Potts chain). 

Our approach is to explore the RFPC by three different 
techniques. We compute the exact GS by using a short- 
est path method^ |, §. This is a generalization of an 
idea that allows to find the GS of the RFIC, and can be 
understood as the zero-temperature limit of the transfer 
matrix computation of the equilibrium spin state. This 



is also employed here, to compare that state with the GS. 
Finally we resort to a qualitative description in terms of 
a random-walk (RW) picture of the groundstate. T his 
is also an extension of the RFIC case, where an earlier 
paper presented a RW description of the exact ground- 
states. 

There is rather little work on the RFPC with the ex- 
ception of trivial field distributions |8|, |j| . Thus it is worth 
to recapitulate some of the features of the RFIC or the 
reasons why it still receives some attention. Due to the 
fact that the model is one-dimensional one could expect 
that the basic magnetization properties are more or less 
trivial. It turns out to be so, however, that in particular 
for binary random field distributions (hi = ±h r ) these 
have very peculiar properties (l0| O). The equation for 
the local magnetization can be considered as a Id dynam- 
ical system with, for suitable parameters, a multifractal 
probability distribution^). For the RFPC this is per- 
haps also the case, generalized to higher dimensions if 
q > 2. In this article we omit such considerations for the 
sake of the equilibrium and groundstate properties. 

One of the points of an exact solution for the GS is 
that it can be investigated how exactly the introduction 
of a small, finite temperature breaks up the GS. For the 
RFPC it is more or less trivial that the overlap between 
the GS and the T > state stays close to one, contrary to 
e.g. spin glasses where this question may still be open to a 
debate. However, it is of interest to study how exactly the 
GS is modified by the "easy" excitations. These are, as 
discussed below, related to the almost degenerate regions 
of the GS. The equilibrium state is also of interest as 
the asymptotic one for out-of-equilibrium processes. In 
such one-dimensional systems the domain walls undergo 
activated dynamics, in particular the RFIC ones perform 
Sinai walks || @. 

The structure of the rest of the paper is as follows. 
First we overview, in section II, the theoretical expecta- 
tions for the groundstate properties and outline qualita- 
tively the consequences for perturbations from a given 
GS, whether by temperature or by changing the local 
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fields (by random shifts or by an uniform applied field). 
Section III is devoted to exact numerical studies of the 
RFPC. In Section IV the temperature is switched on, 
and the effect of entropy on the domain structure is ana- 
lyzed. Section V contains a brief numerical study of the 
effect of correlated fields. Finally, Section VI finishes the 
paper with conclusions. The numerical techniques are 
discussed in the Appendices. 



II. RANDOM WALKS: DECOMPOSITION OF 
THE GROUNDSTATE 

We have earlier presented a way to divide the sequence 
of random fields hi for the RFIC in such a way so that one 
can understand the ensuing domain structure Q. The 
idea is to look at trial random walks: start from a test 
site, and follow the sum of the random fields left /right (an 
arbitrary choice) . Each of such trials constitutes either an 
absorbing or non-absorbing one. The terms describe the 
fact that trials are random walks with absorbing bound- 
aries. One boundary results in an excursion of the RW 
which is such that enough Zeeman energy (2J) is accu- 
mulated to form a domain, but if the other one is met 
first the try is finished and a new one is started. Since 
the RW's are independent, the GS factorizes. Mathemat- 
ically such excursions are given by the following rules. An 
absorbing excursion (ae) is a sequence S of spins start- 
ing at some lattice site i and ending at j > i, with the 
field sum \J2ies ^il for the first time becoming greater 
or equal to 2 J. On the other hand a sequence S' of spins 
from i to j > i is a non-absorbing excursion (nae) if 
< &J2i=i h[ < 2J Vi < k < j ; where a = ±1 is the 
orientation of the spins within the preceding absorbing 
excursion. 

As shown in Figure |l| for each domain one has at least 
one such absorbing RW (it can contain more than one 
such large fluctuation) and it will be bounded by the 
next opposite absorbing RW's on each side. The rules to 
be followed are: (1) Determine an absorbing excursion Sq 
for a given field configuration. If it starts at site io and 
ends at jo, and a is the sign of its field sum, then = a 
V/c G Sq. (2) Starting from jo + 1 find all non-absorbing 
excursions until the next absorbing excursion Si (from i\ 
to ji) is found, whose field sum is by definition opposite 
in sign to the preceding one. The sites k belonging to 
the non-absorbing excursions have the same orientation 
<7fc = (T as those within iSo- The orientation of the spins 
at sites I within iSi is opposite to the later one, o~i = —a. 
(3) Starting at ji + l the search (2) for the next absorbing 
excursion then leads to the overall GS. 

This picture explains also, for binary random fields 
hi = ±h, the degeneracy of the RFIC. For our purposes 
it is more important to note that for any kind of pertur- 
bations it implies that there are three processes, domain 
wall-shifts, destruction and creation of domains, that fol- 
low naturally from the sequence of RW's. If one changes 
the local random fields some of the absorbing excursions 
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FIG. 1: Random walk picture for the RFIC domains. The 
upper part of the figure shows the field sum hi at each of the 
spin sites (starting from the left) and the lower part depicts 
the corresponding domain structure. The random walk starts 
a non — absorbing (nae) excursion as soon as it leaves the 
shaded region of height 2J(for details see text). 

become non-absorbing, and thus domains are destructed 
while new ones may ensue if some non-absorbing ones 
become absorbing. These are in principle easy and in 
practice cumbersome (since one needs to be concerned 
with the exact first passage properties of the perturbed 
field distributions or random walk step length distribu- 
tions) to compute for each process. Moreover, the normal 
outcome is that the domain-walls are shifted. The ma- 
jor results that follow for the RFIC from the RW picture 
concern the domain structure and optimization. The do- 
main magnetic (Zeeman) energy scales linearly with do- 
main size in contrast to the Imry-Ma argument . The 
domain size distribution is exponential, and for /i r < 1 
the Imry-Ma result (Id) ~ l/^r ^ s recovered since both 
the contributions from the non-absorbing walks and the 
absorbing ones that contribute to the domain length scale 
with l/h^ Thus for the RFPC it will be interesting to 
compare both the domain structure and the Zeeman en- 
ergy, 

H z = J2 h i(q$<r-im-l)' (2) 

lev 

For the Potts model one can use the same thought 
experiment: assume the spin a,i = q" , and look for 
"absorbing excursions". Now, in contrast to the q = 2- 
case the general RFPC is more complicated. For each 
of the other q' ^ q" one can follow, analogously to the 
RFIC, the random walk that results from three steps: 
1) null step (local field aligned to other direction than 
q' and q"), 2) positive and 3) negative steps. One has 
q — 1 such trial walks, with the elementary step direc- 
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tions being correlated or shared. These follow "partial" 
sums of local random fields, according to the rules 1) - 3). 
Thus an analogy of the RFIC results in studying the first 
passage properties of a q — 1 dimensional space. Either 
one of the coordinates is increased, or all of them are de- 
creased (g"-step) at the same time. One is interested in 
the first-passage through one of the sides of a hyper-cube 
in q — 1 -dimensions, ie. across the plane where q' = 2 J. 
The walk starts in the immediate vicinity of the corner 
(0, . . . , 0) of the hyper-cube. The simplest complication 
that arises over the q — 2, RFIC case is that for each trial 
walk that compares any of the spin orientations and the 
original one there are "empty" steps consisting of local 
fields aligned with the other q — 2 possibilities. This im- 
plies that the typical domain length might simply scale 
with q, as borne out later by the numerics. 

There are two further complications compared to the 
RFIC, both resulting from the fact that the optimiza- 
tion can not be done as straightforwardly as in the RFIC 
due to the factorization of the landscape to absorbing 
and non-absorbing excursions. Namely for q > 2 the 
process is slightly non-local. Consider the q = 3 case, 
and assume that spin i has q' = 1. If a large fluctuation 
is found by following the appropriate partial sum of the 
q' = 2-orientation, and the next one is of q' = 3-kind, 
one has to consider whether it is energetically more fa- 
vorable to create a 2-domain, a 3-domain, or both, since 
the number of domain walls created can be lowered by 
omitting in this example the 2-domain. One can also use 
the argument the other way: if a large 3-fluctuation ex- 
ists, it is possible to create a 2-domain since one of the 
domain walls is free. The following structure 11113 
3 3 of total energy Ei, is changed into 1112333 of 
energy Ef by flipping a single spin 1^3 (the local ran- 
dom field hi — h directed along a — 3 at the flipped site). 
The latter configuration would be more favourable than 
the former one if, Ef < Ei, yielding h > J. This shows 
that the minimum Zecman energy of a domain of length 
1 should at least be J(q — 1). 



III. GROUND STATE PROPERTIES 

A. Correlations among the successive domains and 
the evolution with increasing disorder 

As noted above the sequence of domains is complicated 
due to the joint optimization. The easiest quantity is 
to compare the probability that the third domain from 
the original one is of the same orientation as the first 
one. We calculate the probability P a (q) of finding the 
every third domain (Pi+2) to be the same or different as 
the first (T>i) one. In absence of any correlations among 
the successive domains one should expect the probability 
P™ c (q) of obtaining T>i + 2 same or different with respect 
to T>i to be, 
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FIG. 2: Left: Probability of obtaining every third domain 
T>i+2 same or of different as the first T>i. The continuous lines 
represents the no-correlation probability function Po C (q). In- 
set : The ratio of P™ c (q)/ P a (q) means the third domain is 
the same as the first one. The system size L = 10 5 and the 
number of configurations 20000. 



( ^ ,for V i = V i+2 
P? C (Q) = { (3) 

We carried out computations with the Shortest Path 
Algorithm (see Appendix A) for q =2, 3, 4, 5 states, with 
the field strength h r = 0.05. It is easy to see that P a (q) 
deviates some from the P™ c (q). Evidently the system 
prefers to have the third domain different from the first 
one. One can put this in two ways: creating a 121 con- 
figuration costs more energy, on one hand, and on the 
other hand as noted before one can add a domain at the 
expense of a single domain wall between two others (131 
becoming 1231). 

The evolution of the domain structure with increasing 
field amplitude h r differs fundamentally from the same 
phenomena in the RFIC. New domains can be created 
in-between two earlier ones at a smaller cost (2 J instead 
of 4 J). Figure H depicts the evolution of the GS domains 
for a q — 5 states sample with identical configuration 
of the random direction of the local fields but increasing 
disorder h r = 0.05, 0.07, 0.15. The large domains are 
split into smaller fragments, mostly at their boundaries. 

B. Distribution of domain length and Zeeman 
energy 

To compare with the q = 2 or RFIC case we now used 
a larger system size of L = 100000, with n = 10 5 samples. 
Figure ^ shows the numerical results for the accumulated 
Zeeman energy Hz as a function of domain length Id for 
different field amplitudes h r . It easy to see that all (Hz) 
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FIG. 4: Log-log plot of the mean Zeeman energy as a function 
of domain length Id- The disorder strength h r — 0.20, 0.15, 
0.10, 0.09, 0.08, 0.07, 0.06, 0.05, 0.04, 0.03 from top to bottom 
in all the subplots. 
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FIG. 3: Evolution of groundstate with increasing disorder 
strength h r for q = 5. 

versus Id curves can be asymptotically (except for the 
smaller Id region) fitted to a straight line of the formQ] 
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FIG. 5: Mean domain length as a function of h r . The straight 
lines are fits to the function f(x) — ax~ 2 where the coefficient 
a =6.99, 6.57, 6.89 and 7.39 respectively for q = 2, 3, 4 and 5. 
The fit works well for smaller h r but deviates for larger ampli- 
tudes. Inset: The average domain length (Id) as a function 
of q. The straight line is a linear fit with slope 201.88. 



(H z ) = qJ + d q h? r l d . (4) 

Here the coefficient d q grows roughly logarithmic for q < 
15. 

Figure |5| shows the mean domain length (Id) for differ- 
ent values of q as a function of the field amplitude h r . The 
data along the x-axis for q = 3, 4, 5 are shifted by fac- 
tors 2, 4 and 8 respectively. In the limit h r <C J the data 
fits quite well with the Imry-Ma scaling|]l5) (Id) ~ l/h%. 
The prefactor of the domain length is (see the inset of 
Figure ||), except for q = 3, linearly dependent on q. Fig- 
ure H finally shows the scaling plots of the probability 
distribution P(ld) of the domain lengths. Apart from the 



initial part, the distribution decays exponentially. We 
may therefore conclude that in spite of the slight correla- 
tions in the domain structure, the RFPC obeys the same 
scalings in the GS as the RFIC. 



C. Response to constant external field 

In the presence of an external field favoring one of the 
spin orientations via an additional term h X^^^i.oi ~ -0 
in the rIamiltonian(dT i arbitrary but fixed and the same 
for all sites), the magnetization becomes of course non- 
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FIG. 6: Scaling plot of the domain length probability distri- 
bution P(ld)- The data are shifted by 7, 14, 21 units along 
the x-axis for q = 3, 4, 5, respectively. Inset: Unsealed plot 
of P(ld) for q = 5. The labels refer to the disorder strength 
h r . 



zero. The ground state magnetization is computed as 

N 

A" 



1 W 



q&ai,a - 1 



1 



(5) 



in a particular direction a (here a £ [1, .., q]) for g differ- 
ent states. Now for a small field strength (h r <C J) in an 
infinite system one has 



rriooih) ~ ft,. 



(6) 



For a finite system the corresponding finite-size scaling 
relation can be obtained as follows |l6|, [17j|. In the GS the 
Zeeman Energy of a block of spins of size L is 



E z ~ L/i. 



(7) 



Further the random field (variance h r ) energy of the 
block is 



Erf L^^h r . 



(8) 



For a non-zero ft one would expect mz,(/i)/i -1 to be a 
function of the dimensionless ratio of Ez and Erf only, 
yielding 



m L (/i) = hm(L 1/2 h/h r ) 



with 



m(x) 



const, for x ^> 1 
1/a; for x <C 1 
Eq. d9|) can once again be written as, 



(9) 



(10) 



rriL 



(ft) = L-WfhiLWh/hr) 



(11) 
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FIG. 7: Ground state magnetisation for q = 2, 3, 4, 5, re- 
spectively. The exponent v — 1.0(for all q) follows from the 
fit to fix) — bx v , which is also verified by the data collapse 
according to Eq. (|lfj|). 



with fh(x) — const, for i»l and m{x) — l/x for iCl. 

The system is subjected to a constant external field 
ft in a fixed direction a — 1 to evaluate m a (S, h) for a 
sample S with some given realization of random fields. 
To start, h = 0.0001, and the GS magnetization is 
computed. Then the field is increased after averag- 
ing 1000 samples, until saturation. Here h r = 0.1, 
L e {5000,10000,20000,40000,80000} while h is in the 
interval [0.0001, 0.03]. For smaller systems one has to 
use higher values of ft r , since otherwise a single domain 
may percolate. 



D. Chaos Exponent 

A slight variation of the random fields in each sam- 
ple gives rise to "chaos". In the case of the RFIC or 
RFIM in general this has already been studied jl^, ^9| . 
There the general behavior is such that for small ampli- 
tudes 8 of the local perturbation the overlap between the 
new and the old GS remains considerable, and the dif- 
ference is a smooth function of 5. This is easy to see if 
one considers the RW factorization of the GS: a variation 
of the local fields gives rise to a new RW, in superposi- 
tion with the original one. In the limit 6 — > this has 
only minute effects except for barely absorbing walks and 
barely non-absorbing ones that can be affected if 8 is large 
enough. Thus the average overlap q is linear for the RFIC 
in 8. Next we look at the chaos properties in the RFPC. 

To modify the initial random fields hi at each site i 
they are replaced by hi = hi + Sh. As before, hi and 
also ft. follow Gaussian distributions and the parameter 6 
measures the strength of the perturbation. Denoting er? 
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FIG. 8: Scaling plot of the overlap correlation function Cs(r) 
versus r8. Data are averaged over 1000 samples for each per- 
turbation strength 8 with fixed system size L = 10 and the 
disorder strength is fixed at h r — 0.05. 



and fff to be the unperturbed and perturbed ground state 
at site i respectively, the ground state overlap correlation 
function up to a distance r is defined by 



C S {r) = 



1 N 

-5>(^ +r )A(of,<?<V r ) 



(12) 



with 



Eh: 




FIG. 9: Destruction of a GS domain for T > in the RFIC. 
The upper part of the figure shows the field sum ]T\ hi in the 
neighbourhood of spin i = 1, 5, the lower part the corre- 
sponding ground state and a T > state. The reason why 
this domain melts at relatively low temperature is indicated 
in the middle part: The field sum of spin 5 is just outside the 
shaded area(see text). 



A(cr J ff / )=2^ 1 « r '-l. (13) 
In the limit N — ► oo one would expect a scaling form 

Cs{r) ~ c(r/m) (14) 

with 

{(fjKt 1 , (15) 

In Figure || we show the results of our calculation of the 
GS overlap correlation function Cs(r). r is varied from 4 
up to 160 for each perturbation strength S. The obtained 
data collapse agrees quite well with the argument above. 

IV. DOMAINS AT FINITE TEMPERATURE 

A. Changes with the introduction of a finite 
temperature 

For T > we calculated the expectation values (<7j) 
using the transfer matrix technique pi B2j as described 



in Appendix B. The domains start "melting", i.e. the 
overlap with the GS, {<Ji<Ji^Gs)i decreases, as the tem- 
perature is switched on from T = 0. The effect of tem- 
perature is similar to a random field perturbation in the 
sense that both should play a role in particular where the 
fluctuations in the random landscape (absorbing/almost 
absorbing walks) make it the easiest. On the other hand 
the effect of temperature is not limited to those regions, 
only. Next we compare the finite temperature and GS 
configurations to illuminate the T > physics. 

In Figure ^| we represent the destruction of a domain 
for a system of size of L= 400 and h r = 0.7, as the 
temperature is raised to 0.2 from T = 0. The random 
field distribution is uniform in the interval [—1,1]. The 
zero temperature domain from i = 2 to i = 5 is due to a 
single absorbing excursion, depicted by the large shaded 
region in the upper part of the Figure |9|. The inset in the 
middle part shows the sensitive region at i = 5, where 
the field sum just crosses the absorbing boundary. As 
the temperature is increased from zero, the entropy gives 
weight to states where the whole region or much of it is 
flipped, to the a = — 1 state instead of +1. 

Figure [To| shows the melting of a GS domain at T = 0.1 
for a system of L = 400 and strength of disorder h r = 1.0 
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T = 0.7. Between T = 0.25 and T = 0.7 AT = 0.05. The 
system size is restricted to L — 400 and h r = 0.5 



FIG. 10: Creation of a fluctuation for T > (RFIC). The 
upper part of the figure shows the field sum Yli hi in the 
neighbourhood of spin i = 1, 10, the lower part the corre- 
sponding ground state and two T > states. The reason why 
this domain melts at relatively low temperature is indicated 
in the middle part: The field sum of spin 4 and 10 are very 
close to the absorbing boundary (see text). 



(the distribution of local fields is again uniform in [-1,1] 
). Now the sum-of- fields random walk in the GS from i 
= 1 to 10 is actually part of a non-absorbing excursion 
inside the embedding domain with a = — 1. It is however 
susceptible to the introduction of temperature (note the 
spins i = 4 and 10 where the RW is closest to creating a 
domain, in the inset). 

At a finite temperature the sites 2, 3, 4 form a "new 
domain" with a = — 1. This new domain contains only 
one part of the RF fluctuation, from 1 to 4. 

To further investigate the melting of a GS domains at 
finite temperatures, an overview of the GS and equilib- 
rium configurations at finite temperatures is shown in 
Figure JO]. The calculations for each temperature in- 
cluding the GS are carried out for the same local field 
configuration. The thick dotted line corresponds to the 
GS whereas the continuous lines represents the melted 
configurations with the gradual increase in temperature. 

In general, for very low temperatures there are only a 
few segments of spins, which differ from their respective 
GS orientations H , growing with increasing temperature 
in height and width. At high temperatures (T> J and 
T 3> h r ) the expectation value (a) of each spin finally 
fluctuates around the mean value l/q. 



B. Domains at T > 

To construct a domain structure at any finite temper- 
ature we focus to the melting of individual q states. It is 



easy to calculate the contribution for each of q states to 
the resultant melting by taking a single state in the diag- 
onal matrix, in the transfer matrix calculation, and keep- 
ing the others zero. The maximum of these contributions 
will give rise to a definite local spin state (ctj = 1, ....q). 
Figure |l~2|(a) represents the GS (straight bold line) and 
melted (<Tj) at finite temperatures for q — 3. Figures 
jl2](b) , (c) , (d) shows the components (for each q) of prob- 
ability P{pi — 1,2,3) at T > of individual spins (<7j) 
at sites i. Consider the 1st segment of domain in (a). It 
is evident that the melted domain approaches the state 
q = 3 as confirmed by the behavior of the same segment 
in Figures (b), (c) and (d). To analyze the melting pro- 
cess quantitatively we investigate the probability distri- 
butions of lengths of the melted segments and the melting 
rates where m, = (<7j) for each of q =2, 3, 4, 5. 

We also compute the finite temperature distribution of 
domain lengths. The simulations are for different tem- 
peratures with a system of fixed L = 300 and disorder 
strength h r — 0.7. The distribution of local random 
fields is uniform in [-1,1]. The disorder, temperature, 
and system size are restricted by the fact that one has 
to compute products of L(q x q) transfer matrices. One 
also prefers to avoid Gaussian disorder instead of one 
with a bounded support. It is clear from Figure O, that 
except for the initial part, the probability distribution 
P{It) varies exponentially with It, similarly to GS do- 
mains. With an increasing temperature the decay in the 
tail becomes faster, indicating that correlations diminish 
as expected. 

We define the probability distribution of melting rates 

Amj 
AT 



as, 



AT 



nL 



with 



Am = \a{T + AT) - a(T)\, 



(16) 



(17) 
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FIG. 12: The GS of RFPC with g = 3 and equilibrium states 
((<Ti) is the thermal expectation value of the spin value at 
site i) for T > 0. Within the interval 0.1 < T < 0.45 the 
temperature is raised by AT = 0.05. (b), (c), and (d) shows 
components of melting probabilities for q = 1 . . . 3. 
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FIG. 13: Probability distribution of finite T domains P(It) 
for q = 2, 3, 4, 5. The data are averaged over 1065 samples. 
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FIG. 14: Probability distribution F(^f-) at constant field 
h r = 0.7. 



where n is the number of samples. We measure the 
change in magnetisation for each spin cr, when the 
temperature changes by AT. For counting we use the 
resolution 6 = 0.001, for which N(^ip) is actually the 



number of spins a l with e - § , 



Am 
AT 



Figure |lj shows the distribution P(^S) at different 
temperatures for q —2, 3, 4, 5; AT = 0.05 is the same 
for all temperatures and 10 5 configurations. The data 
points for are not shown as in this regime we observe 
a strong maximum at — > 0. For q — 2 there is a small 
but distinct peak in the curve for T = 0.2 which gradually 
diminishes for higher values of q. 
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FIG. 15: Scaling plot of the average T > domain length for 
different h r - The data are shifted along the y-axis by -0.1 and 
-0.2 units for q =3, 4 respectively. L — 300 and the data are 
averaged over 50000 configurations. 



C. Finite temperature lengthscales 

Finally we investigate how the local regions of the mag- 
netization evolve at finite T@. We calculate the average 
length l m (l m — (It)) from the finite temperature con- 
figurations. Figure |l5| represents how this lengthscale 
changes with temperature, if one first scales away the 
T = dependence on the field. Note that the temper- 
ature is restricted so that T > 0.09. A further collapse 
with the right combination of h and T makes it possible 
to observe an universal scaling function for l m so that 



h v 

l m = ld,F( — ) 



where the scaling function F — > 
nent 

v 2/3 



(18) 

1 as T -» 0. The expo- 
(19) 



does not change with q, and thus one has the same scaling 
for q > 2 as for the RFIC. The implication is basically 
the same: one has per each domain a number of "easy" 
excitations that allow the magnetization to change a lot 
from the GS. The presence of such almost degenerate 
regions is not affect by the exact value of q. 



V. RFIC WITH CORRELATED RANDOM 
FIELDS 

We have already seen that for an uncorrelated random 
field distribution the average domain length asymptoti- 
cally follows the Imry-Ma argument, i.e. (l^) ~ l/h~ 2 . 
This should change for spatially correlated random fields. 
Now we obtain an exponent 7 with (Id) ~ Vv 7 where 
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FIG. 16: Scaling exponent 7 of p(ld) = h' 1 P(KJ.ld) against p. 
p is the exponent of the correlated random field distribution 
(see text). 



7 approaches 2 as the spatial correlations decrease. Con- 
sider a power-law correlation function of the form 



(hih 



i+r / 



(20) 



To use an analogy of the Imry-Ma argument, the random 
field energy of a domain of length Id is expressed as 



Erf oc 



\ 



hih i 



(21) 



Making use of the correlation function (21) we readily 
obtain, 



E 



RF 



\ 



hi 



l-p 



(22) 



with < p < 1. The domain wall energy is as usual 

Edw — J (23) 

To create domains, in the GS, Erf and Eow are of the 
same order of magnitude, thus 



Id — K 7 



where the exponent 



7 



(24) 



(25) 



for p < 1. For p > 1 this Imry-Ma type argument predicts 
the irrelevance of the correlations in the disorder, i.e. 
p = 2. 

We computed the probability distribution P(ld) on a 
large system of size L — 65536, with the local, correlated, 
Gaussian random fields generated by the Noise Construc- 
tion Algorithm pfj|. Statistical averaging was performed 
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over 10 5 samples. As in the case of uncorrelated fields 
P(ld) is again exponential for all values of p we studied. 
The disorder h r was varied for each value of p, to calcu- 
late the exponent 7 from the data collapse of P{ld) ver- 
sus la- 7 deviates significantly from the predicted value 
as p approaches unity, but finally approach the value 2 
for p — > 00. The deviation might be due to logarithmic 
corrections to the Imry-Ma prediction ( |25| ) close to the 
critical value p = 2. The same properties can also be ver- 
ified for the RFPC, if the random field magnitudes and 
their directions obey the same correlations. 



VI. CONCLUSIONS 

It is a fortunate fact that one can solve numerically, but 
exactly, for the groundstates of the RFPC. For the RFIM 
one has access to powerful graph optimization methods 
that work in d arbitrary, but for q > 2 it can actually 
be shown that the problem is NP-complete for d > 1 
pl| . Here we have used this to study the Potts chains, 
augmented with random walk considerations and with 
transfer matrix calculations. 

As it is already known that the RFIC groundstate is, 
essentially, separable to independent regions and what 
the consequences are it is most natural to compare the 
q > 2 behavior to the RFIC physics. It turns out that 
essentially all more general scaling behaviors, whether 
pertaining to the GS or to the finite temperature states, 
follow similar laws. This is true, in groundstates, for the 
Zeeman energy, for the form of the domain distribution, 
and for the average domain length. At T > we have 
demonstrated that the change from the GS configurations 
follows a similar course for q arbitrary. Finally also for 
GS chaos the value of q has no essential importance. 

There are of course slight differences since the local do- 
main structure has correlations due to the "decorations" 
of large domains with smaller ones that utilize the fact 
that only one domain wall needs to be created if e.g. in- 
duced by increasing the disorder strength. These are of 
secondary importance for the scaling properties. 

The greatest deviations from the expected behavior we 
observed in the case of the RFIC with correlated fields, 
where the scaling exponent of the domain length with 
disorder strength does not seem to take the expected 
value 7 = 2 (but is slightly less) for weakly correlated 
fields. This may be due to corrections to scaling, or to 
the detailed properties of a RW picture, similarly to the 
uncorrelated case, when the walks are fractional Brown- 
ian motions. 

In the case of the RFPC it is possible to propose a 
number of topics to study. Eg. the generalizations of 
the multifractal magnetization distributions of the RFIC 
seem mathematically interesting. Also the detailed prop- 
erties of the "melting" of the GS might merit further 
study. 
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FIG. 17: The graph corresponding to the Shortest Path 
Method for the RFPC problem. 
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APPENDIX A 

The ground state of g-State RFPC can be calculated 
exactly by mapping onto a shortest path problem on an 
acyclic graph. The graph G(L, q) corresponding to a q- 
State RF Potts chain of length L consists of (L+l)xq 
nodes which we consider to be arranged like presented in 
figure 1 and denoted by (x, k) with x = 1, L and k = 
1, q. Each node (a;, k) is connected to all of its right 
neighbours, i.e. to all nodes (x + 1, k') with k' = 1, ....,(7 
by a directed one-way edge with cost 

C( x ,x+i : k : k>) = -J(l S k,k' ~ X ) - h x+i(l S k,& ~ !)• ( A1 ) 

Fixing the spin at x = to be in state k the groundstate 
of RFPC then corresponds to a shortest path from (1, k) 
to (L,k) in G(L,k). To apply periodic boundary condi- 
tions it is necessary to solve q shortest path problems, 
one for each of the q (with (q, q) spins at sites x = 1 and 
x = L+l) possible states. One of these shortest paths 
with minimum energy gives rise to the optimal ground- 
state. 

Since the underlying graph (Fig. [l?]) is acyclic no cy- 
cle, in particular none with negative weight, can occur 
and Dijkstra's algorithm for finding the shortest path in 
a weighted graph can be applied although not all edge 
weights are positive. The running time of this algorithm 
is linear both in number n of the nodes and in the num- 
ber m of the edges, thus the overall complexity to find 
the groundstate of the RFPC with periodic boundary 
condition is 0(q 2 L). 
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APPENDIX B 



It is easy to see from Eq. ([l7]) that 



For temperature T > we generalize the Transfer Ma- 
trix Method H| for the g-States RFPC. The calculation 
of the partition function can be reduced to the prob- 
lem of finding the product of N random matrices: 



JV 



Z N = J2 exp{-/3Ha) = Tr JJ L« 

{a} Vi=l / 



(Bl) 



where 

r(i) . 



[ exp[/3J(g^- !<Tl - 1) + (3hi{q5„- jt<Si - 1)] ] 



(B2) 



Here {j,k = l,...,q). 

Using this expression of partition function the expec- 
tation value (a r ) for each spin can be expressed as, 



(a r ) 



(B3) 



(oy) 



Tr 



n f =lL (i + -i) S 



(B4) 



with Sjj = z.<5y. 



The only computational effort consists in calculating 
the product of N q x q transfer matrices. Since the ele- 
ments of Lw become very small for low temperatures, ar- 
bitrarily small temperatures can not be considered. How- 
ever, the admissible temperature interval is sufficient for 
our investigations. 
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